Hybrid Charmonium and the p — tx Puzzle 
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Using the method of QCD sum rules, we estimate the energy of the lowest hybrid charmonium 
state, and find it to be at the energy of the ^'(25') state, about 600 MeV above the J/*I'(1S') state. 
Since our solution is not consistent with a pure hybrid at this energy, we conclude that the ^''(25') 
state is probably an admixed cc and hybrid ccg state. From this conjecture we find a possible 
explanation of the famous p — tt puzzle. 
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I. INTRODUCTION 

A hybrid meson is a state composed of a quark and 
antiquark color octet, along with a valence gluon, 
giving a color zero particle. Hybrids are of great 
interest in studying the nature of QCD. The non- 
perturbative method of QCD sum rules [T] has long 
been used to predict the energies of light quark hy- 
brid mesons [2| and hybrid baryons (see Ref. [sj). 
In the present work we use this method to estimate 
the energy of the lowest charmonium hybrid. A ma- 
jor motivation for the present work is to understand 
the nature of the ^''(25') state, and to find a possible 
explanation of the long-standing p — tt puzzle. 

The p — TT puzzle concerns the branching ratios for 
hadronic decays of the '^'{2S) state compared to the 
J/'^{1S) state. By taking ratios of hadronic decays 
to gamma decays for these heavy quark states the 
wave functions at the origin cancel, and by using 
the lowest-order diagrams one obtains the ratios of 
branching rates for two charmonium states 



R = 



B{'f'{2S) -> h) 



B(J/*(15) 
0.13 , 



h) B{J/^{1S) -> e^ 



(1) 



the so-called 13 % rule. For the '^'{2S) state com- 
pared to the J/^'(1S') state, however, the hadronic 
{e.g.p — tt) decay ratio is more than an order of mag- 
nitude smaller than predicted Q • This is the p — n 
puzzle. There have been many, many theoretical at- 
tempts to explain this puzzle: Chen and Braatan Q 
review earlier work by Hon and Soni, Brodsky, Lep- 
age and Tuan, Karl and Roberts, Chaichian and 
Tornqvist, Pinsky, Brodsky and Karliner, and Li, 



Bugg and Zou. All seem to agree that this is an un- 
solved puzzle. More recently there has been an at- 
tempt to locate the source of the problem with 
the suggestion that there is a cancellation of two 
processes in the ^'(25*) decay. Our present work 
suggests that one can obtain such a cancellation by 
including valence gluonic structure. 

In Sec. II we show that the energy of the lowest 
hybrid charmonium state with the quantum num- 
bers J^^ = 1 is that of the *'(2S'), but our so- 
lution is not consistent with a pure hybrid. Since 
the cc(2S') state is also expected to have that energy 
(about 600 MeV above that of the J/^'(1S') state), 
we predict that the 5*' (25*) state is an admixture of 
cc(2iS') and hybrid components. In Sec. HI we show 
that this can provide a solution to the p — tt puzzle. 
In Sec. IV we give our conclusions and compare our 
results to lattice gauge calculations. 



II. HYBRID 1" CHARMONIUM USING 
QCD SUM RULES 

We now will use the method of QCD sum rules 
to attempt to find the lowest hybrid charmonium 
state, assuming that such a pure hybrid charmonium 
meson with quantum numbers J^'~^ — 1 exists. 
First, let us review the method and the criteria for 
determining if one has obtained a satisfactory and 
accurate solution. 



A. Method of QCD Sum Rules for a Hybrid 
Charmonium Meson 

The starting point of the method of QCD sum 
rules is the correlator, which for a hybrid meson is 



u^'^ix) = (r[j^(x)jj^(o)]) 



(2) 



with the current J'^{x) creating the hybrid state be- 
ing studied. The QCD sum rule is obtained by eval- 
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uating n'"' in two ways. First, after a Fourier trans- 
form to momentum space, a dispersion relation gives 
the left-hand side (Ihs) of the sum rule: 



/Ihs 



g2) 



, Imn^^fs) 

ds- y- 3 

7r(s — g^j 



where Ma is the mass of the state A (assuming zero 
width) and Sq is the start of the continuum-a param- 
eter to be determined. The imaginary part of n(s), 
with the term for the state we are seeking shown as 
a pole (corresponding to a 5{s — M\) term in Imll) 
and the higher-lying states produced by Jjj shown 
as the continuum, is illustrated in the figure: 



Im[n(s)] 




fit to continuuin 



FIG. 1: QCD sum rule study of a state A with mass Ma 
(no width) 



Next n''^(g) is evaluated by an operator product 
expansion (O.P.E.), giving the right-hand side (rhs) 
of the sum rule 



(4) 



with increasing k corresponding to increasing dimen- 
sion of Ok- 

After a Borel transform, B, defined in Appendix 
B, in which the q variable is replaced by the Borel 
mass, Ms, the final QCD sum rule has the form 



1 lAi^ „ f°° , Ini[n(s)l 

I" Jsa 7r(s-g^) 



Bj2^k{q){0\Ok\0) 



(5) 



as we shall show below. Note that the Borel trans- 
form produces an exponential decrease with increas- 
ing values of s, as shown by the pole term in Eq®. 
This reduces the contribution of the continuum, 
where s > M^. 

This sum rule is used to estimate the heavy hy- 
brid mass, Mhhm- One of the main sources of er- 
ror is the treatment of the continuum. In addition 



to the parameter So, one must parameterize the ef- 
fective shape of the continuum. The criteria for a 
satisfactory solution are: 1) The contribution of the 
continuum should not be as large as the pole term 
in the Ihs. 2) With an exact sum rule, the value of 
Mhhm is independent of the value of Mb- With 
the approximation of a fit to the continuum, there 
should be a minimum or maximum in the value of 
Mhhm vs Mb, and the value of Mhhm at this ex- 
tremum should be approximately the value of Mb- 
3) There should be a gap between the solution for 
^^HHM ^'^'^ ■^o' which reduces the contribution of 
the continuum to be smaller than the pole term due 
to the factor of e~'^I^B after the Borel transform, as 
we shall explain below. If the value of So is much 
larger than the expected excited hybrid states, how- 
ever, the solution is not physical. 



B. Previous Results for Hybrids Using QCD 
Sum Rules and Lattice QCD Methods 

Since mesons with certain quantum numbers, such 
as 1 ^, cannot have a standard qq meson composi- 
tion, there is a strong motivation for both experi- 
mental and theoretical searches for hybrid mesons, 
which can have such states. These states are called 
exotic or hermaphrodite mesons. Shortly after the 
introduction of QCD sum rules, they were used to at- 
tempt to find the masses and widths of exotics. The 
most accurate calculations [1| predicted 1 hybrid 
light-quark mesons in the 1.3-1.7 GeV region, where 
such a state has been found . The solutions satisfy 
the criteria for a good solution, as explained above. 
For example, for a solution that has a hybrid mass 
of 1.3 GeV, the value of So was 1.7 GeV, which is a 
reasonable separation of the lowest from the higher 
hybrid states. On the other hand, for these light me- 
son hybrids many terms in the O.P.E. are needed, 
which significantly increases the uncertainty. This 
explains why previous calculations have found a 
rather wide range of values for the light-quark exotic 
hybrid. 

For heavy-quark hybrids the higher-order terms 
in the O.P.E. are quite small, so the QCD sum rule 
method is more accurate. Since we are not trying 
to predict exotic hybrids, however, there is the seri- 
ous complication of nonhybrid meson-hybrid meson 
mixing, which will be discussed in detail below. This 
is an important aspect of our present work. 

There have been many lattice QCD calculations 
of glueballs and light-quark hybrids [1]. The most 
recent calculations of light-quark hybrids find the 
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lightest exotics to be about 2 GeV [9j, quite a bit 
higher in energy than the QCD sum rule calcula- 
tions. This probably is due to the fact that lattice 
QCD calculations for light quarks have some incon- 
sistencies at the present time, while they are much 
more accurate for heavy quark systems [l0| , as are 
QCD sum rules. Exotic charmonium states have 
been calculated using lattice QCD, and the 1 ^ was 
found to be about 4.4 GeV 11;], with the expecta- 
tion that the 1 hybrid charmonium state is at a 
similar energy. For the calculation of nonexotic hy- 
brid mesons, such as 1 hybrids, there are other 
difficulties for both methods, as we shall discuss be- 
low. 



C. Heavy Hybrid Meson Correlator 

For a hybrid meson with quantum numbers 1" 
we use the standard current 2] 



(6) 



with — C71/, where C is the charge conjuga- 
tion operator, 7^ is the usual Dirac matrix, and 
\E' is the heavy quark field. Carrying out a four- 
dimensional Fourier transform, the correlator in mo- 
mentum space is 



n^''(p) 



(27r)4 
Tr[G^^G''^)]{pi) , 



(7) 



where S""^ is the quark propagator, with colors a 
and b. The color properties of C?^", the gluon 
color field, with /.t, a Dirac indices, are given in Ap- 
pendix A. Note that the traces are both fermion 
and color traces. Details oiTr[S''^TaS^''T p\{p) and 
Tr[G^''^G''l^)\{pi) are given in Appendix A. 

It is important to recognize that the correlator 
used in the QCD sum rule method is similar to the 
correlator used in the lattice gauge approach, with 
the same objective of finding the mass of a heavy 
hybrid meson. 

As described in the preceding subsection, the cor- 
relator is evaluated in the method of QCD sum rules 
via an O.P.E. and a dispersion relation. 



D. O.P.E. of the Scalar Correlator in 
Momentum Space 

The QCD sum rule method uses an operator prod- 
uct expansion (O.P.E.) in dimension (or inverse mo- 



mentum), Eq(|3]). For the hybrid meson the lowest- 
order diagram jl2l | is shown in Fig. 2: 




FIG. 2: Lowest-order term in sum rule 



Using the standard quark and gluon propagators 
(see Appendix A for some details) we find for this 
lowest-dimensional process 

nr(p) = -6 / 

^ Pi 



. (27r)4 
Pi - (j>- Pi 



-2 WYi+P^iP-PiY (8) 

Pi 



-2p>t{p-pir]}{ 



2 M, 



Q 



i{p-piY '4 



+ [-AP''Pif 



with 



Jo a(l 



da 



a)p^ - Ml 



(9) 



Extracting the scalar correlator H'^, defined by 

n^'^(p) = ip^p./p" - gn^^{p) + {p^.Pu/p^n^{p). 

and carrying out the pi integrals, one finds for 
nf (p), the scalar term of li'^'^ [p), 



nf(p) 



da 



-{ 



d(3p^p 



i^^rJo a(l-«)7o p^(l-/3)-^ 

4 155 Ml 319 , 

3a- a^ 12 (a-a2)2 + 12^ Q,_a2 

iMy-|M|/]+ (10) 



4 /•! 



-1 



2 -/o "'■'p2(i_;3)__^ 3 

^){-3P + pys) - (4/3 



a ~ a 

^2 I a3 



3/3^-^/373)]} 
terms with three and four integrals . 
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We find that the terms with three and four integrals 
are very small, and we do not include them in our 
calculation. 

The second term in the O.P.E for the heavy-quark 
hybrid correlator includes the gluon condensate, il- 
lustrated in Fig. 3: 




FIG. 3: Gluon condensate term in sum rule 



For this process, the correlator Ht^'^ (p) has the 
same form as Eq. ([4]), except that the gluon trace 
used to obtain H^"' {p) is replaced with the trace over 



IS repi, 

the gluon condensate [l|, which gives 

-g^^gn ■ (11) 

The fermionic factor is the same as for Fig. 2. From 
this one finds 

2 Anr2\T , ^M^i^r 11„2 53 2 



[p' - 4A/^)/o + p'Yi^p' - -jMq 



4 4 

6 p^ A p'^ 

giving for the scalar part of Fig. 3 
3i < > 11 



(12) 



n^(p) 



2 {An\ 



-|m|p^ + Hm^)/o. (13) 



E. Borel Transform of the Correlator 



with higher-order terms being very small. We shall 
see that even 112 is essentially negligible within the 
accuracy of the method, and the convergence after 
the Borel transform is established. 

Using the equations in Appendix B, from 
Eqs. I|37l37l38p . where the quantities Bi through Bg 
are given, with BIlf{p) = flf (M^), we find 



fif {Mb) = -^^[-WM^B,-155Ml^B2 

+319B3 - IQM^Bi - 164^5 
-110(^6 - M^Bj) - IGM^Sg] (14) 



1 



-(1+5) 



{[-310^^ + 638(5 - 656(1 -f 5)] 
i^3(2||(l + ^)) + [-1892^4 
606(5 - 1968(1 + (5) - 32(4(5 - 3 



l + (5 



-^)]^2(2^(1 + .)) 



-4778- 



1 + 5 



Ml 

-2 

9442(5- 4920(1 -^(5) 



-128(4^-3 



6^ 



6^ 



1 + 5 3(1 + 5)- 



r)] 



i^i(2^(l + ^)) 

-t-f-1356 



(15) 



-96(45-3 



M2 



1 + 5 

(52 



6284(5 - 3280(1 -I- (5) 

(53 

r)] 



l + ,5 3(1 + (5)- 



ii'o(2-— ^(1 + (5)) + multiple integral terms . 

The multiple integral terms in Eq.(12) are small and 
are dropped, and 5 is a variable of integration. 

In a similar way, taking the Borel transform of 
nf (p), with Snf (p) EE nf (Mb), we find 



To ensure convergence of the O.P.E. of the corre- 
lator, one performs a Borel transform B 1], defined 
in Appendix B. This is discussed in detail in the 
early papers on QCD sum rulesfT], For our prob- 
lem we have assumed that 11"^ (p) ~ nf (p) -I- Ilf (p), 



tlUMs) 



2(47r) 



(16) 



14 M2 M2 
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F. QCD Sum Rule for Hybrid Charm Meson 

The method of QCD sum rules uses a dispersion 
relation for the correlator, which it equates to the 
correlator's operator product. Following the usual 
convention we call the dispersion relation the left- 
hand side and the operator product expansion the 
right-hand side: Ilihs = dispersion relation, H^hs — 
operator product expansion. Neglecting the width 
of the hybrid meson, the dispersion relation has the 
form of a pole and a continuum: ds . The 
dispersion relation is evaluated in Euclidean space, 
—Q^, and the continuum is assumed to start 
at s = So- After the Borel transform, the form we 
use for the Ihs is 

(17) 



xe 



with F the numerator of the pole, and Li and L2 
constants used to fit the form of the continuum. We 
shall use a standard method to eliminate F, and fit 
the sum rule requirement that the solution should 
not be sensitive to Mb ■ In theory, the exact solution 
should be independent of AIb- 

For convenience in carrying out the sum rule, we 
have fit the rhs of the correlator to a polynomial in 
the Borel mass, IIi{Mb) = aiAf| -t-aaM^ -f asA/l + 
a4M| and fl2(Ms) = 6iAf| 4- feaM^ + hMl^ + 
h^M% with ai = 122.29, aa = -143.88, as = 
45.79, a4 = -0.346, 61 = 616.0, 62 = -279.1, 63 = 
-226.56, and 64 = 144.515, with units GeV^. 

We find that, for all values of Mb relevant to the 
sum rule, 112 is just about 1 % of 111, and drop it, as 
the method is only valid to a few percent. The sum 
rule is obtained by taking the ratio of 11 to diij^j-ili. 
To do this we use the relations 



-T- A 



v>0 



-Ae 



B 



2" 



A 



A 



A 



A 



d,,Mie "'^Koi-j^) = -Ae [ifo(-^) 

B B 

-AV-^)l. (18) 

Taking the ratio of ^l/]^J2^Illhs{MB) = 
c^i/MS nrhs(MB) to the equation nihs(MB) = 



Ilrhs{MB), one obtains the sum rule for the mass of 
the heavy charmonium hybrid meson: 



M 



HH = {e '"'b[so{LiMI + L2M%) + LiM% 

+2L2M|]-fai/M|nf} (19) 



L2M*: 



The result of the QCD sum rule fit is shown in 
Fig. 4, with so=60.0 GeV2, Li=-99.0 GeV\ and 
^2=6.06 GeV^: 



> 

D 

13.7 



CM 



13.6 
13.5 
13.4 
13.3 



13.1 13.2 13.3 13.4 13.5 13.6 13.7 
(m/ (GeV)2 

FIG. 4: QCD sum rule for heavy hybrid charm meson 

The parameters Li and L2 have magnitudes ex- 
pected for the fit to the continunmm. Except for 
the large value of So, the solution satisfies the crite- 
ria for a successful QCD sum rule: as Fig. 4 shows, 
the solution for the mass is within a few per cent of 
the value of Mb in the region of stability. The mass 
predicted for the hybrid is 



Mhhm ^ 3.66GeV 



(20) 



within a 10% accuracy of the sum rule method, while 
the experimental mass of the ^''(25') state IIS] is 



M{-9'{2S)) = 3.68GeV 



(21) 



The large value of So, however, predicts that the 
excited hybrids are at a very high energy, as found 
in lattice gauge calculations, and that the ^''(25') 
cannot be a pure hybrid. Therefore, we expect that 
the physical ^'{2S) is an admixture with a charm 
meson and a hybrid charm component. A second, 
orthogonal mixed state will be in the continuum. 
As we now show, this can provide a solution to the 
p — TT puzzle. 
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III. HYBRID-NORMAL CHARMONIUM 
AND A POSSIBLE SOLUTION TO THE p - -k 
PUZZLE 

In our treatment of hadronic decays of the ^''(25') 
state, we use the Sigma/GluebaU model, which was 
motivated by the BES analysis of glueball decay [31 
and the study of scalar mesons and scalar glueballs 
using QCD sum rules [M, lb\ . We briefly review this 
model below. 




FIG. 5: Scalar glueball-meson coupling theorem 



A. The Sigma/Glueball Model 

In energy regions where there are both scalar 
mesons and scalar glueballs, it is expected that 0'''+ 
states will be an admixture of mesons and glueballs. 
For this reason, when using QCD sum rules to find 
such states one must use currents that are a lin- 
ear combination of glueball and meson currents. A 
scalar glueball current can have the form 

J<=(x) = asG'' , (22) 

while a scalar meson current has the form 

J"(2:) = ]^{u{x)u{x) - d{x)d{x)) . (23) 

We use for our 0++ current [3] (with Mo needed for 
correct dimensions) 

Jo++ = /3M„J,„ + (1-1/3 I) Jg. (24) 

The QCD sum rule calculation makes use of the 
correlator n(x) = {T[Jq++Jq++]). The cross term 
between J™ and J*^ is evaluated by using the scalar 
glueball-meson coupling theorem Il6|: 

j dx{T[J^{x),r^{Q)]) ^ < > , (25) 

with < qq >= the quark condensate , which is illus- 
trated in Fig. 5: 

The results of the QCD sum rule calculations [ij, 
[Tsj l are that there are three solutions: 

80% scalar glueball at 1500 MeV the fo(1500) 

80% scalar meson at 1350 MeV -> the fo(1370) 

Light Scalar Glueball 400-600 MeV ^ the 
Sigma/Glueball 

The Sigma/Glueball model follows from the exis- 
tence of the a, a scalar vr — tt resonance with a broad 



width and the same mass as the scalar glueball found 
in the sum rule calculations, and makes use of the 
glueball-meson coupling shown in Fig. 5. It has 
been used for the study of the Roper resonance de- 
cay into a nucleon and a cr [1, [13] , and the prediction 
of a production in proton-proton high energy colli- 
sions [18], and other applications. The relevance for 
our present work is that, just as it is possible to de- 
termine the mixing parameter for a state consisting 
of a scalar meson and a scalar glueball, it should be 
possible to determine the cc and hybrid ccg admix- 
ture for charmonium systems, and for T systems. 

B. Hybrid Mixing Model and the p — n Puzzle 

The solution to the p — n puzzle by the mixing of 
hybrid and normal meson components of the ^'(25) 
state can be understood from Fig. 6 and Fig. 7, 
which illustrate the decay of a cc and a ccg state 
into two hadrons. The cc decay involves the matrix 
element < 7rp|0|'I''(cc, 25*) >: 




FIG. 6: Lowest-order PQCD diagram for a cc decay into 
two hadrons 



The corresponding hybrid decay involves the ma- 
trix element < TTp\0'\'if'{ccg,2S) >, with the dia- 
gram shown in Fig. 7. 

Assuming that the 2s state is a cc-ccg admixture 

\^'{2S)) - b\'f'{cc,2S)) + Vl-b^\'f'{ccg,2S)) , 
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FIG. 7: Decay of a ccG state into two hadrons 



and recognizing that the O and O' matrix elements 
are approximately equal, we see that the solution to 
the p — TT puzzle requires 

b+^/l~b^ ~ 0.1 , 
The solution of the p — n puzzle would be given if 
b ~ -.7 . 

In other words, if 6 ~ —.7, we have found a solution 
of the p — TT puzzle. The evaluation of b is rather 
complicated, and the testing of this conjecture will 
be carried out in future work. 



IV. CONCLUSION 

Using the method of QCD sum rules, we have 
shown that the ^''(25') state cannot be a pure char- 
monium hybrid. We have found that the energy of 
the lowest J^'-' = 1 hybrid charmonium state is 
approximately the same as the ^''(25') state, about 
600 MeV above the J/^'(1S') state, but that the 
QCD sum rule solution is not consistent with a pure 
hybrid. The standard model prediction for cc{2S) 
is at approximately the same energy. Therefore we 
expect that the physical ^''(25') state is an admix- 
ture of a cc{2S) and a 00(8)5(8) (25). Using this pic- 
ture,we find a possible solution to the famous p — n 
puzzle. 

There have been many lattice calculations of ex- 
otic hybrid mesons. There is experimental evidence 
for an exotic light-quark 1 meson (see Ref Q for 
a discussion) at 1.4 to 1.6 GeV, which is consistent 
with QCD sum rule calculations 0, while lattice 
calculations find the lowest 1 hybrid at 1.9 to 
2.1 GeV. The lowest energy 1 ^ charmonium hy- 
brid found in lattice calculations is at 4.4 GeV (see 
Ref [11]), about 800 MeV above our 1~ hybrid 
charmonium solution. This is consistent with our 
large value of Sq- For the 1 state, we have shown 



that one must use a mixed cc and ccg current, which 
we shall use in future work. This use of a mixed cur- 
rent to define the correlator has not been done in lat- 
tice QCD caluclations, but with the QCD sum rule 
method it can be done in a rather straight-forward 
calculation. 

It is interesting that the energy difference between 
the T(2iS') and the T(IS') is also approximately 600 
MeV. If this is the energy of a T{nS) hybrid, this 
could provide a solution to the puzzling a decays of 
T{nS) states that have recently been observed 
Investigation of this system is a topic of future re- 
search. 



V. APPENDIX 

A. Heavy-Quark Hybrid Correlator in 
Momentum Space 

The current to create a heavy-quark hybrid meson 
with J^"^ = 1 is 

= ^T^G'"'^ , (26) 

where ^ is the heavy quark field, T„ = Cj^, 7^ is 
the usual Dirac matrix, C is the charge conjugation 
operator, and the gluon color field is 

G"" = Ev^'"' (27) 

Q=l ^ 

with Xa the SU(3) generator {Tr[XaXb] — 26ab)- 
From this one finds for the correlator of a heavy- 
quark hybrid meson 



n'^-(p) = 




''^■^ <T[J^ix)J^m> (28) 



Tr[G''^G''^)]ipi) , (29) 

where is the quark propagator, a standard Dirac 
propagator for a fermion with mass Mq , with colors 
a and b. Using G = 47^7'', TrljaJts] — 4,gaf3, and 

Tr[7A 7a 757/3] = '^{gXagSp " g\S9al3 + g\l39Sa)- 

{p\ks ~ k\ks){g\agsi3 - gxsga/3 + gxpgsa)] 

[k^-Mmp-kY-Mi] ■ (^°) 
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To use dimensional regularization we define the 
quantity D = 4 — e, and let e ^ to complete the 
integrals. One can then show 



Two key equations which we need are (with Ki, 
the modified Bessel functions) 



1 



(27r)^ ik^-Mmp-ky-Ml] 



with 



lo = 



da 



{a - a^)p^ - Ml 



,(31) 



(32) 



(fc-l)!(Af|)'=-i 



The first term in the operator product expansion, 
shown in Fig. 2, has a standard gluon propagator. 
For the gluon trace (see, e.g., Ref [3]) one finds 



^,_i^l_fc, ^ 2{yy/^K,{2V^) . (37) 
b 



Transforms used in the body of the paper are 
^ da 



Tr[G^"G'^'^](pi) = -47r^z(5a/3^ 

Pi 



Mi/ a 

PlPl _L PlPl 



Pi 



PfPl 

-5m/3 2~ 

Pi 



(33) 

Pi 



The correlator for the process shown in Fig. 2 is 
found to be 



(27r)- 
Pi ■ ip-Pi) 



P\ 



[p^'Pi+P^iP-PiT 



9 A/^ 

-2p'^ip-p,r]}{- " 



-.2 + hip ~ P^^'^ 



41 2 



1 



r)]/o(p-Pi)}(34) 



6 Q 3 Q(p_p^): 

The next term in the operator product expansion 
for this heavy quark system, where quark conden- 
sates are negligible, is the gluon condensate term, 
shown in Fig. 3. The trace over the quark propaga- 
tors is the same as in Eq. (l30l) . The gluon field trace 
for this term is 

< Tr[G^'^G-'0]{p,) > = {27t)^S\p,)^ < > 

i^gt^-g'^P - g'-^gn • (35) 
From this one finds 112"', given in Eq. (12). 



B. Borel Transforms 

A key method that enables one to use the operator 
expansion to get accurate sum rules is the use of the 
Borel transform i3, defined by 



B 



lim _-L_(q2)"(_^- 
q^,n^oo (n- 1)!^^ ' dq^- 



■ (36) 



Ba 



-B{p'-Ml/{a-a')y 



(a-a2)' 
= 2e"2*^«/*'-f«ii:o(2A/A/Af|) 



Bi ^ B 



da 



d(3P 



[a-a^Y 7o n-/?^fn2 ^Ik 1 

U PhP (a-Q2)(i_0)J 



2Af;^ 



8 -2 



Af 



/f2 



db—.e [2X2(2^(1 + J)) 



1 + 5 



B 



Af2 A'/2 
+8X1(2^(1 + 5)) + 6ifo(2^(l + 5))] 

B B 



B 



= 2Af^ 



da 



' S -2^(1+5) Af^ 



d(3l3 



1 + 5 

V /I I f\\ I on i^^ /o V / 



Af 



-f 12X2(2^(1 + 5)) + 30K,{2j^{l + S)) 
Af 2 

+20fCo(2^(l + 5))] 



B 



da 



df3P 



2Af^ / dSSe 
Jq 

Af2 



-2— S-(l+<5) Af^ 

[2K,{2-^{1 + S)) 



12f^2(2^(l + 5)) + 30f^i(2^(l + 5))] 



Ml 



+20K^{2j^{l + 5))] 



dpp 



1 

(Q-a2)(l-,3)J 



2Af^ / d55e 
Jo 



_2^(l+<5) Af^ 

Af2 



[2f.2(2^(l-f^)) 



Sr. = B 



da 



(--^ PnP {a-a^)(l-l3)\ 



9 



2M6 / d55{l + 5)e ^ [2X^(2-^ {1 + 6)) 

Jo ^^'^B 



MA M% Bs = B ^p' 



+ 12K2{2-^{1 + S))+30K,{2-^{1 + S))] Jo Jo {1 - p)[p 



da 4 d/3(4/3 - 3/?2 + /53/3) 



52 53 -2^(1+.-) +8ifi(2-|(l + <5)) + 6ifo(2-|(l + <5))].(38) 
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